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Impedance matched thin metamaterials 
make metals absorbing 
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^ AEgisTech., Nanogenesis Division 41 0 Jan Davis Dr, hHuntsville, AL 35806, USA, "^Dept. of the Army, Charles M. Bowden Facility, 
Redstone Arsenal, AL 35898, USA. 

Metals are generally considered good reflectors over the entire electromagnetic spectrum up to their plasma 
frequency. Here we demonstrate an approach to tailor their absorbing characteristics based on the effective 
metamaterial properties of thin, periodic metallo-dielectric multilayers by exploiting a broadband, 
inherently non-resonant, surface impedance matching mechanism. Based on this mechanism, we design, 
fabricate and test omnidirectional, thin (<1 micron), polarization independent, extremely efficient 
absorbers (in principle being capable to reach A > 99%) over a frequency range spanning from the UV to the 
JR. Our approach opens new venues to design cost effective materials for many applications such as 
thermo-photovoltaic energy conversion devices, light harvesting for solar cells, flat panel display, infrared 
detectors, stray light reduction, stealth and others. 



Absorbing materials are crucial in many applications such as thermo-photovoltaic energy conversion 
devices^ light harvesting for solar cells^, flat panel displays^ stealth, and stray light reduction in imaging 
devices^. In nature, black pigments are found in the compound eyes of most insects to reduce crosstalk^. 
Another case is the black color on the Papilio Ulysses butterfly wings that is due to a combination of black 
pigment and structured surface that traps the light^. Many venues have been proposed during recent years to 
achieve broadband, polarization -insensitive, omni- directional absorbers. Currently available broadband low- 
reflectance materials include gold black^, black nickels, including porous nickel phosphorus (NiP)^ and black 
paint. Black metals and paints tend to have rough, protruding surface features to trap the light causing it to 
undergo several reflections. A material with 15% reflectance can absorb >99% of the energy after just three 
reflections. Black materials with rough surfaces such as anodized aluminum perform best at near normal 
incidence and degrade significantly for large angles^. Recently, successful demonstration of an extremely dark 
material has been made from low density nano-tube arrays^"^\ However, fabrication of mechanically and 
environmentally durable nano-tube coatings over large areas is still a challenge and involves costly processes. 
Emerging new technologies have focused on absorption mechanisms based on the excitation of localized and/or 
surface plasmons on nano-patterned metal surfaces. For example, in Ref. 12 the authors experimentally dem- 
onstrate an ultrathin plasmonic absorber over the entire visible spectrum, in Ref. 13 the authors numerically 
study a broadband, 1-D metal grating with almost perfect absorption at }i ~ 600 nm, in Ref. 14 broadband 
absorption is predicted for all-metallic 1-D and 2-D gratings, in Ref. 15 an absorber based on light propagation 
in a metamaterial forming an effective "black-hole" is discussed, in Ref. 16 it is proposed a negative-index based, 
wide angle absorber for IR radiation, and in Ref. 17 broadband absorbers and selective emitters based on 
plasmonic Brewster metasurfaces are numerically investigated. 

In this work we propose a viable, yet effective, alternative to achieve efficient absorbers spanning from the UV 
to the near IR by matching the impedance of a metamaterial to the impedance of air. The geometry of the 
multilayer stack under consideration is depicted in Fig. 1. It consists of an AT-period, symmetric structure whose 
elementary cell is made of a metallic layer sandwiched between two identical dielectric layers. The metamaterial 
does not require surface features or nanopatterning and can be deposited over large areas with traditional 
deposition techniques. 

In Ref. 18 we have shown that in general any stratification with a symmetric elementary cell but not limited to a 
metallo-dielectric stratification, or a symmetric three-layer elementary cell described in Fig. 1, is equivalent to an 
effective bulk metamaterial, having the same physical thickness of the equivalent multilayer and univocally 
described by an effective characteristic impedance (Z^^) and the Bloch wave vector (Kg) of the elementary cell. 
The expression of Z^fp Kb and the free space characteristic impedance (ZyJ are respectively: 
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Figure 1 | Generic metal/dielectric/metal symmetric elementary cell of 
period A = Id^ + dM- The overall structure consists of N repetition of the 
elementary cell grown on a substrate. 
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where A is the period of the symmetric elementary cell, oo the fre- 
quency, niij the elements of the transfer matrix of the elementary cell, 
Zo the vacuum impedance, and 0 the incident angle. Additional 
details of the derivation of Eqs.(l) and (2) are presented in the 
Methods. We note that our homogenization procedure, different 
from conventional approaches, does not rely on the approximation 
that the incident wavelength must be much longer than the dimen- 
sion of the elementary cell. It maintains a general validity for any 
periodic stratification, not necessarily metallo -dielectric, character- 
ized by a symmetric elementary celP^. The propagation length (Lp), 
i.e. the distance for the intensity of an electromagnetic wave prop- 
agating inside the metamaterial to decay by a factor of can be 
extracted from the imaginary part of the Bloch vector: Lp = 
(2Im(KB))~\ The propagation length (Lp) together with the effective 
characteristic impedance (Z^jj) are the key quantities to understand 
the absorption mechanism in these structures. Ideally, in order to 
achieve an efficient absorber we need a metamaterial that satisfies 
two requirements: first, it has to possess an effective characteristic 
impedance that matches the free space one for all incident angle and 
all polarization (Zg^ = Zfs) and, second, it has to possess a very short 
propagation length, say few elementary cells. The first condition, Z^^ 
= Zfs, ensures that all the electromagnetic energy impinging the 
surface of the metamaterial enters into it with no reflections, while 
the second condition ensures that all the energy is efficiently 
absorbed within a few elementary cells. Bulk metals could in prin- 
ciple be good absorbers as they exhibit a consistent imaginary part of 
the permittivity, but actually they are not because, below the plasma 
frequency, the electromagnetic radiation practically cannot pen- 
etrate into them due to their high impedance mismatch, and thus 
they reflect most of the light. In what follows, using metamaterial 
concepts, we demonstrate both theoretically and experimentally a 
method to tailor their effective impedance and match it to the 
vacuum over a broad frequency range which can span the Vis and 
near IR. 



For a metallo -dielectric stratification with the symmetric three- 
layer elementary cell described in Fig. 1, the effective impedance at 
normal incidence from Eq.(l) and (2) can be explicitly written as: 



(3) 



/ sin (dMkonM)sD + Inunp tan (np/co^p) cos (^m^o^^m) - tan^ (nDkodp) sin (^m^o^^m) 
y sin (dMkonM)sM + ^UMno tan (nDkodo) cos (dukoriM) - &d tan^ (nDkodo) sin (dMkonM)'' 

where nM and no are respectively the refractive index of the metal 
and dielectric layer, while Sm and Sd their respective relative electric 
permittivities. Here for simplicity we analyze theoretically the nor- 
mal incidence case, nevertheless all the resulting considerations can 
be straightforwardly extended to any arbitrary incident conditions 
and polarizations. For thin metallic layers, i.e. kodM< <1, a first 
order expansion of the trigonometric functions in (3) gives the fol- 
lowing expression for the effective impedance: 



1 dukoSM tan^ (riDkodD) - 2nD tan (^d/co^^d) - dukoSn 



Zo rij^ y dukoSD tan^ (riDkodD) - Iud tan {nDhdo) - duhsM 
Eq.(4) can be recast in a more convenient form as: 
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where the numerator function f^M and the denominator function 
are two parabolas respectively given by: 



fN{x)=X^SMS —IUdX — SdS, 
fD{x)=X^£DS -IHdX-SmS, 
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with x= tdLn{nDkodD) and 3 = kodM< <l- For the sake of simpUcity, 
in the analysis that follows we neglect the imaginary part of the 
permittivity of the metal (as we will see later on, the results that stem 
from this analysis hold also when realistic metallic dispersions are 
used). Obviously, when the radicand of Eq.(5.a) is negative the effec- 
tive impedance will be a pure imaginary number and therefore the 
metamaterial will show metallic-like behavior with 100% reflectivity 
and complete impedance mismatch, while in the case that the radi- 
cand of Eq.(5.a) is positive the effective impedance will be positive 
and the structure will show a dielectric -like behavior. In order to 
study the sign of the radicand, Eq.(5.a) can be rewritten, using 
Ruffini's theorem, as: 



8m (x-X-n) (x-X+n) 
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where the zeros of/jv(:>c) and/o(x) are respectively given by: 
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For 3>l/ ^J — Em (Case I) Ufe is rather dull, in fact the radicand of 
Eq.(6.a) will always be negative and the structure will show a metal- 
lic-like behavior. After all, this is the regime one would expect given 
the presence of metal layers. Luckily, things become much more 
interesting for ^ < 1/ ^ — Sm (Case II). In this latter case, the radicand 
becomes positive when either x < x < x+m or X-p<x<x+p and 
in those intervals the structure will show a dielectric-like behavior In 
Fig. 2 we show a graphical representation that summarizes the 
results. 
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Figure 2 | Graphical representation of the study of the sign of the radicand of Eq.(5.a). In Case II bands of metallic-like behavior alternate with bands of 
dielectric-like behavior. 



In particular, Case II of Fig. 2 shows the bands where the structure 
either behaves as a metal or as a dielectric. The two dielectric regions 
are bounded respectively by the zeroes and the poles of the effective 
impedance. Namely, the dielectric region X-n <x<x+n is boun- 
ded by Zeff^Zo = 0 (zeroes of fN(x)), while the dielectric region 
X-D<x<x+D is bounded by Z^^Zq = oo (zeros of/o(x)). A further 
analysis of the effective dispersion of the conditions Z^ff/Zo = ^ and 
Z^yZo = 0 shows that they can be approximately expressed as: 



Zo 

Zpff 



-- 00 





Ism , , rn 

— flM+ ^ 

2 8d 2nD 


m 


= 0,1,2, 


\^dD = 


^^^+^ 
2 AfiD 


m 


= 0,1,2, 


dD= — 


1 8m J 2m - 1 

2 8d ^flD 


m 


= 1,2,3, 


dD= — 


1 m 

- + — ^ 

2 2nD 


m 


= 1,2,3, 



(7) 



dielectric-like behavior, but, in principle, it is also possible to 
match its impedance to the vacuum by a proper choice of the 
dielectric material, no matter how many metal layers the structure 
is made of. The structure could contain hundreds of metal layers 
and still be impedance matched with vacuum! 

In Fig. 3 we plot, as an example, Sd vs. X/X^ according to (9) for 
metal layers with different thickness whose dispersion we model 
by using a generic, free-electron gas, lossless Drude model 

^M = '^ — {^/^pY i^p = 27i;c/(j0p is the metal plasma wavelength). 

In Fig. 3, it can be seen that the values of the permittivity of the 
dielectric materials necessary to achieve perfect impedance matching 
are within the range of the permittivity of many dielectric and/or 
semiconductor materials. For example, in Fig. 4 we show the ref- 
lectance at normal incidence i^= 1 1 — (Zg^/Zo)/l + (Zg^/Zo) | 
numerically calculated in the plane (d^/^^p, X/Xp) for djvi/^p = 1/30 
and Sd = 2.25. As expected from the previous study, the structure 



Let us now analyze in particular the dielectric region X-D<x<x+Dy 
i.e. the one bounded by Z^^Zq = oo. The impedance at the center of 
this band is: 
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By solving Eq.(8) under the condition of perfect impedance matching 
with vacuum, we find: 



'Z^ 
Zo 



■■io8d= — - — ^ — -0 — -• (9) 

1 + 8mS^-S^ ^ ^ 



Eq.(9) dictates what is the best choice of the permittivity of the 
dielectric material for perfect impedance matching in the dielec- 
tric band. Therefore, not only our structure possesses bands of 




Figure 3 | £d vs. X/X^ according to (9) for a lossless free electron gas. 
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Figure 4 | Reflectance in the (do/^p, for djvi/^p = 1/30 and £d = 

2.25. The metal layer dispersion is described with a generic, free-electron 
gas, lossless Drude model. 

shows two bands of dielectric-like behavior. The one that is bounded 
by the poles of the effective impedance shows perfect impedance 
matching. The white dashed lines represent the dispersion of the 
poles and the zeroes of the effective impedance calculated according 
to (7). 

Now that we have analyzed in detail the properties of the effective 
impedance of our metamaterial and the impedance matching con- 
dition, we come to the analysis of the propagation length (Lp). As 
already mentioned in the introduction, a good absorber should be 
characterized by a short propagation length. At normal incidence for 
a semi-infmite (s-i) structure the absorption is given by the well-known- 
formula: As-i = ARt{Zeff/Z^) I (^\Zeff/Zo\' + 2Re{Zeff/Zo) + iy By 

expanding the Bloch vector at the first order for thin metal layers, 
it can be demonstrated by simple mathematical manipulations that, 
at the center of the pass band of the Bloch vector, i.e. for cos(KbA) = 
0, Lp can be approximated as: 



A lm(8M)kodM ' 

Eq.(lO) tells us that the propagation length is inversely proportional 
to the imaginary part of the metal permittivity (dissipation), as one 
may expect, and directly proportional to the refractive index of the 
dielectric material. In order to have an efficient absorber a metal with 
high dissipation properties, such as Cu for example, would be desir- 
able in conjunction with a low refractive index dielectric material, 
ideally as close as possible to vacuum. So far we have concentrated 
for the sake of generality on metals with nonspecific dispersion 
properties. We want now provide two examples in which metals 
with realistic dispersion properties are used. In particular, we con- 
centrate on copper (Cu) and silver (Ag). In Fig. 5, we plot the bulk 
dispersion relations of the two metals, as reported in the book of 
Palik^^, it is worthwhile to note that copper shows intrinsic better 
absorption properties, i.e. higher imaginary part of the permittivity, 
with respect to Ag. In order to be as realistic as possible, we have also 
taken into account the increased resistivity of the thin metallic layers 
when their thickness becomes comparable or shorter than the elec- 
tron mean free path through the Fuchs-Sondhemeir theory^°'^\ More 
details about the procedure used to evaluate the dispersion of the 
thin metallic layers can be found in the Methods. 

In Fig. 6 and 7 we show the numerically calculated absorption 
at normal incidence for a semi-infinite structure (A^./) and Lp/A 
in the (X.,dD) plane, respectively for a nD(dD)/Cu(10 nm)/nD(dD) 




A. (|Lim) 

Figure 5 | Bulk dispersions for Ag and Cu as reported in the book of 
Palik^^ 

elementary cell and for a nD(dD)/Ag(10 nm)/nD(dD) with n^ = 
1.52 (Si02). In both Fig. 6(a) and Fig. 7(a) broad (non-resonant) 
regions of impedance matching and high absorption (A > 99%) 
are evident. Those regions are approximately bounded by the 
poles of the effective impedance (white dashed lines), consistent 
with the results of the previous study summarized in Fig. 2. Bands 
of dielectric-like behavior are also present and they are approxi- 
mately bounded by the zeroes of the effective impedance, as 
expected. Note that, although in this case we fully take into 
account the metal dissipation without neglecting the imaginary 
part of its permittivity, nevertheless the general picture previously 
described still holds, as expected. 

Fig. 6(b) and 7(b) show Lp/A in the Ck^di)) plane. The yellow line 
represents here the center of the pass band calculated according to 
the condition cos(KbA) = 0. It is seen that Lp in the Ag case is much 
longer than the Cu case. 

In view of Eq.lO, the reason for the different performances is 
ascribed to the intrinsic better absorption properties of Cu, i.e. higher 
imaginary part of the permittivity, with respect to Ag, as it is seen 
from the bulk dispersion properties of both materials reported in 
Fig. 5. Fig. 6(c) and 7(c) compare Lp/A calculated at the center of 
the pass band according to the exact formula Lp = {2Im(KB))~^ and 
according to (10). 

The results on the propagation length suggest that the choice of 
the particular kind of metal is essential to obtain efficient absorption 
in short distances. While in fact, in the examples made, both struc- 
tures ensure a high degree of impedance matching, the Cu-based 
structure is much more efficient in absorbing the e.m. radiation than 
the Ag-based structure. This means that, based on the respective 
propagation lengths, a 6 -period Cu-based structure would be highly 
absorbing whereas the same structure made of Ag would be highly 
transparent. Although here we are interested in increasing the 
absorption properties of the structure, it is nevertheless interesting 
to note the great versatility of the impedance matching mechanism 
that we have elucidated: it makes possible to pass from a highly 
absorbing structure to a highly transparent structure just by a proper 
choice of the particular kind of metal. 

Results 

Based on the theoretical study presented in the previous section, we 
have fabricated two samples: sample 1 and sample2 (see Methods for 
the fabrication details). In both cases the elementary cell is a three 
layer stratification with 10 nm of Cu sandwiched between two layers 
of Si02 30 nm thick. Sample 1 is made of 6 repetitions of the ele- 
mentary cell for a total of L ~ 3Lp, sample2 is made of 10 repetitions 
of the elementary cell for a total length L ~ 5Lp. Both samples have 
been grown on a quartz substrate. The transmittance (T) of the two 
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Figure 6 | (a) Numerically calculated absorption of the semi-infinite structure in the (Xydj)) plane for a nD(dD)/Cu( 10 nm)/nD (do) elementary cell with 
Ud = 1.52 (Si02). The white dashed lines represent the poles and the zeroes of the effective impedance calculated according to (7). Bands of impedance 
matching and almost unitary absorption bounded by the poles of the effective impedance are visible, as predicted, (b) Lp/ A in the (Xydj;)) plane. The yellow 
line represents here the center of the pass band calculated according to the condition cos(KbA) = 0. (c) Lp/A vs. X calculated at the center of the pass band 
according to the exact formula Lp = {2Im{KB))~^ and according to (10). 



samples has been measured between 0° and 70° over a wide spectral 
range (300 nm up to 800 nm) while the reflectance (R) has been 
measured between 20° and 70° over the same spectral range (see 
Methods for measurements details). In Fig. 8 and Fig. 9 we report 
(left panel) the measured average transmittance (T) and reflectance 
(R) (the average has been taken over TE and TM polarization) as 
well as the derived absorption (A = 1-R-T) for samplel and sample2 
respectively in the (}t,0) plane, in comparison with the same 



quantities as predicted by the theory (right panel). A close look at 
Fig. 8 and Fig. 9 reveals that the reflectance from the two samples is 
almost identical, while the transmittance is slightly higher for sample 
1 and scales approximately as T ~ exp(—L/Lp) in agreement with the 
theoretical predictions based on the propagation length analysis. As 
expected, both samples show high degree of polarization independ- 
ent absorption (>90%), over a wide range of incident angles and a 
spectral band extending from UV to near IR. Note also that the first 



SCIENTIFIC REPORTS | 3 : 3203 | DOI: 1 0.1 038/srep03203 



5 



[SiO2(dD)/Ag(10nm)/SiO2(do)] 




0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 

X (|Lim) 




0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 

A. (|Lim) 




0.75 



X (|Lim) 

Figure 7 | Same as in Fig. 6 for the nD(dD)/Ag(10 nm)/nD (do) elementary cell. 



sample has a total length of only 420 nm and a total metal thickness 
of just 60 nm. The experimental absorption near 750 nm is larger 
than predicted by theory and we attribute this to the possibility that 
the 10 nm thick Cu film is below the percolation threshold which 
could induce losses due to the excitation of localized surface 
plasmons. 

Discussion 

We have reported a practical and simple way to achieve broadband, 
omnidirectional, polarization-insensitive absorption using periodic, 
thin, metallo- dielectric, multilayers in their metamaterial operation 
regime. Those structures can be easily fabricated over large areas at 



low cost using standard sputtering or thermal evaporation tech- 
niques^^. We remark that, although in this work we have focused 
on Si02/Cu multilayers for absorbing purposes, different combina- 
tion of metal/dielectric materials can be exploited. As an example, in 
Fig. 10 we present the theoretical results for an N = 10 -period struc- 
ture whose elementary cell is made of a tungsten (W) layer 5 nm 
thick sandwiched between two Ta205 layers 75 nm thick. In this case 
the structure shows an ultra- wide band of abosrption ranging from 
300 nm to 2.5 jim. 

The possible applications for these structures are innumerable. 
Absorbing materials structured at the nanoscale level are a fun- 
damental part in the design of any efficient thermo-photovoltaic 
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[SiO2(30nm)/Cu(10nm)/SiO2(30nm)]^6 




Figure 8 | (Left Panel) Average transmittance, reflectance and absorption measured for sample 1 in the plane. (Right panel) Average transmittance, 
reflectance and absorption numerically calculated. 



energy conversion device or solar cefl. We also remark that the 
mechanism exploited in this work is based on a non-resonant imped- 
ance matching phenomenon which is therefore independent of the 
particular nature of the substrate. This makes possible, for example, 
to use our thin metamaterial to cover bulky metaflic structures to 
eliminate electromagnetic reflection over an extremely broad fre- 
quency range from UV to almost mid-IR. As an example, we have 
grown sample 2 on a reflecting surface made of an optically thick 
silver layer. A picture of the sample and measured absorption is 
shown in Fig. 11. 

Moreover, the non-resonant nature of the phenomenon also allows 
us to use these metamaterials in conjunction with other known 
absorbing configurations, such the classical Salisbury-mirror-like 



configuration, i.e. antireflection coating/absorber/mirror, to further 
enhance their performance. 

We also mention that the same approach that here we have 
employed to obtain efficient absorbers metamaterials for the visible 
and near-IR could be employed as well to obtain efficient absorbers 
in the mid-IR by using polar material such as SiC, for example. 

Methods 

Metamaterial description of multilayers. The fields at the input and output surface 
of a single layer material (a Fabry-Perot etalon) are linked to each other through a 
matrix relation that contains information on the angle of incidence by the 
longitudinal momentum of the impinging wave, the material properties 8 and |i, 
and the layer thickness d^^ (see Fig. 12): 
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Figure 9 | Same as in Fig. 8 for sample 2. 
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where di is the thickness of the i* layer, k^ is the longitudinal component of the k 
vector in the i* layer, a-^ = i^j for TE polarization or oCi = Si for TM polarization, being 
|ii and Si the bulk magnetic permeability and the electric permittivity of the i*^ layer. 
For a symmetric structure^*^ it is possible to define kf ^ and oCeff so that: 



For TE polarization, a = ^ = E, (p = (l/|i)dE/dz. For TM polarization, oc = 8, ^ = 
H and cp = (l/s)dH/dz. In a cascade of layers, the fields at the input and output 
interface are still related each other through a matrix relation. In this case the transfer 
matrix of the whole structure is the product of all the single layer transfer matrix. If we 
consider an elementary cell of length A = ^di we have: 



/ Cos(kf A) 



^S/.(kfA)\ 



(13) 



(POut 



n 

Elemetary Cell 



( Cos(^di) ^Sm(k^di)\ 



^ --iSm(k^di) Cos(k;di) 



(12) 



Note that in general Eq. (13) provides 4 equations for 2 unknown quantities and in 
general there is no solution, unless the elementary cell is symmetric. 

In Eq. (13), oCeff = i^eff (effective permeability) for TE polarization, or oCeff = Sgff 
(effective permittivity) for TM polarization, kf ^ is the longitudinal component of the 

k vector in the effective medium: kf ^ = Kb = neff (w/c)^ — (kx/ngff)^ where, kx is 
the transverse component of the k vector, that is continuous through the interfaces, 
k^ is also equal to the Bloch vector of the structure^*^. 
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Figure 10 | Absorption calculated for a 10 period Ta205(75 nm)/ 
W(5 nm)/Ta205(75 nm) structure. 

From Eq. (13) we can derive 

^Cos(kf A) = m^\ = m^2 



The characteristic impedance in a generic medium takes the following expressions^"^ 
Z' = \ l- T- f or TE polarization 



£ kz 
' s k 



(15) 



for TM polarization 



Eq.(l) follows from Eq.(14) and (15). Finally, for N repetition of the elementary cell 
the transfer matrix is given by: 



V 



Cos(kf iVA) 



-^SinikfNA) 



Cos(kfiVA) 



(16) 



In other words, the N-period structure is equivalent to a bulk of the same effective 
medium of the elementary cell and having the same length L = NA of the N-period 
structure itself: i.e. the dispersion properties of the single, symmetric elementary cell 
and the dispersion of the N-period structure made of the same elementary cell are 
identical. 

Dispersion properties of thin metallic layers. Metal layers with dimensions on the 
order of the electron mean-free path (/) or less, exhibit higher resistivity than bulk 
conductors. The effect was first observed experimentally by Thomson in 1901^^ and 
later theoretically explained by Fuchs in 1938^° and Sondheimer in 1952^\ According 
to the Fuchs- Sondheimer theory, the thin layer conductivity {(^thin) is related to the 
bulk conductivity (a), through a function of the layer thickness (d) over the electron 
mean-free path (^ = d//): 



(y thin 



O(0 



with 
1 



12 ■ 



^'•(-0-^(1- 



5 1 



where Ei( — Q is a special function known as exponential integral: 

fco -t 

^'(-{) = -J 



(17) 



(18) 



(19) 



For a Drude metal, with plasma frequency cOp and damping y, the electric permittivity 



' takes the following expression: 



co{(o-\-iy) 



(20) 



1 



The damping term is inversely proportional to the conductivity y oc -, therefore 

o 

we can derive a thin layer damping (ythin) in analogy with the thin layer 
conductivity: 

ythin ^ ^ ^ ^(0 
y (Tthin ^ 



(21) 



i^'^) and a thin layer Drude dispersion for the electric permittivity e^^" 



Drude _ 
^thin - 



coi (o + iy 



(22) 



Unfortunately, close to the plasma frequency most metal, including Cu and W, do not 
strictly follow a Drude dispersion. In order to take into account the increasing 
absorption of thin layers, we subtract from the experimental bulk data taken from the 
book of Palik^'^ ) the free-electron contribution i^^^fif^) and we subsequently 
reintroduce the free-electron contribution of the thin layer (ej^^jj'^^) calculated 
according to Eq.(22): 





0.5 0.6 

X([im) 

Figure 1 1 | [Left] Measured absorption of a Si02/Cu multilayer stack (sample 2) on a reflecting surface made of silver and [Right] picture of the sample 
taken outdoors under normal daylight. The reflection of the blue sky and clouds is apparent in the picture on the reflecting side of the sample. 
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= n 

ElemelaiyCell 



di d 

' A 

^Sinik\d^) Cos(k;d,) 



Figure 12 | Transfer matrix representation for a single layer and a generic three-layer structure. 



Hulk 



Drude _ 
■ Hulk 



(23) 



Note that in the limit of thick layers, i.e. d ;» /, Eq.(23) gives the correct limit to the 
experimentally measured data of the bulk dispersion. The electron mean free path for 
Cu is 39 nm and for W 19 nm at room temperature, according to data reported in 
literature^^'^^. 

Fabrication. The metal-dielectric stacks were fabricated using a Denton magnetron 
sputtering system. The Si02 film was deposited using a 3" Si02 target and the copper 
layer was deposited using a 3" copper target. A 700 micron thick quartz wafer was 
used as a substrate. The system can be used to sputter a 6" diameter substrate with 
uniforminity of about 5% across the substrate. 

Optical measurements. Spectral transmission and reflection measurements were 
done using a WooUam VASE spectroscopic variable angle ellipsometer. The system 
allows measurements for s (TE) and p (TM) polarization from 15 degrees to grazing 
angle of incidence in reflection and 0 to grazing in transmission. Transmission 
measurements were taken with respect to no sample in place. Angular measurements 
were done from 20 to 70 degrees with an increment of 5 degrees. 
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